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In this paper we investigate the universality and scaling properties of the well-known quantities 
in classical statistical mechanics near the quantum phase transition point. We show that transverse 
susceptibility and derivatives of correlation functions with respect to the parameter that drives 
the quantum phase transitions, exhibit logarithmic divergence and finite size scaling the same as 
entanglement. In other words the non-analytic and finite size scaling behaviors of entanglement is 
not its intrinsic properties and inherit from the non-analytic and scaling behaviors of correlation 
functions and surveying at least the nearest neighbor correlation functions could specify the scaling 
and divergence properties of entanglement. However we show that the correlation functions could 
capture the quantum critical point without pre-assumed order parameters even for the cases where 
the two-body entanglement is absent. 
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I. INTRODUCTION 

Quantum Phase Transition (QPT) is a continuous 
phase transition which occurs at zero temperature where, 
quantum fluctuations play the dominant role. Suppres- 
sion of the thermal fluctuations at zero temperature in- 
troduces the ground state as the representative of the 
system. Often the abrupt changes in the ground state 
manifest itself in the non-analytic behavior of macro- 
scopic properties of the system. But in a number of 
important cases this criterion fails to predict a QPT, 
such as the three-dimensional metalinsulator transition 
of non-interacting electrons in a random potential. The 
QPT are typically accompanied by certain exponents 
pertaining to the class of quantum phase transitions 2 . 
The concept of correlation, i.e. information of one part 
of a system about another part, is a key element in 
(QPT) physics. Correlations can emerge from classical 
and quantum sources. 

The existence of genuinely quantum correlations can 
be usually inferred by the presence of entanglement 
among parts of a system. Indeed, entanglement dis- 
plays a rather interesting behavior at QPT, being able 
to indicate a Quantum Critical Point (QCP) through 
nonanalyticities^— . Moreover, for one-dimensional criti- 
cal systems, ground state entanglement entropy exhibits 
a universal logarithmic scaling^. In several years an idea 
has arranged that a complete classification of the critical 
many-body state requires the introduction of concepts 
from quantum information theory^. 

Although it is known that the ground state and corre- 
lation functions could capture the QCP, but to the best 
of our knowledge, universality and scaling properties of 
well-known correlation functions have not been investi- 
gated. In this work, we will study the universality and 
scaling properties of correlation functions and the Trans- 
verse Susceptibility (TS) to compare with universality 



and scaling properties of entanglement. We show that 
derivative of correlation functions with respect to the 
transverse field, exhibit a logarithmic singularity at the 
QCP and show the scaling behavior close to the QCP the 
same as TS. The scaling behavior and constants which 
are obtained in this work are the same as that obtained 
in refs. [3[ and [J] for the entanglement. In other words, 
divergence of entanglement close the QCP is not surpris- 
ing since the reduced density matrix and consequently 
the two-body entanglement are smooth functions of cor- 
relation functions and its non-analytic behavior and finite 
size scaling rise out of divergence and scaling behavior of 
correlation functions' derivative. However, we show that 
the correlation functions between far neighbors are able 
to characterize a QPT, even for distances where pairwise 
entanglement is absent. This is a consequence of the 
longer range of correlation functions at the critical point 
in comparison with the short-range behavior of pairwise 
entanglement. 



II. XY MODEL IN TRANSVERSE FIELD AND 
NEAREST-NEIGHBOR CORRELATION 
FUNCTIONS 

The system under consideration is the antiferromag- 
netic spin-1/2 XY chain in a transverse magnetic field, a 
model central both to condensed-matter and information 
theory and subject to intense studyi^. 

The hamiltonian is 



N 



H = jy x afat+ 1 + J y aVoV +1 + hal 



(1) 



where J x > and J y => are the exchange cou- 
plings, h is the transverse field and N is the number of 
sites. This model is exactly solvable^ with the Jordan- 
Wigner transformation which in the momentum space 
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leads to H = J2k 2e k(ltlk - 5), where 7^(7*;) is the 
fermion creation (annihilation) operator. For states of 
even fermions, 



£fc = V [h + ( J x + J y )cos(fc)] 2 + (J a - J y ) 2 sm\k), (2) 



and the ground state energy is given by Eq = 
n = - 1. For J„ 



with fc 
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equa- 
tion (JTJ reduces to the Ising model in transverse field 
(ITF), whereas for h = it is the XY model. For all 
the interval < J x + J y < 1 the models belong to 
the Ising universality class and for N — 00 they un- 
dergo a quantum phase transition at the critical value 
A c = (Jx + Jy)/h = 1. The magnetization < a x > is 
different from zero for A > 1 and it vanishes at the tran- 
sition point. On the contrary the transverse magnetiza- 
tion < a z > is different from zero for any value of A. At 
the phase transition the correlation length £ diverges as 
£~ |A C - A|- y , with v = lii. 

The two point nearest neighbor (NN) correlation 
functions and transverse magnetization are coincidently 
the expectation value of the coupling term in the 
Hamiltoniani^, 



G' 



dJ x 



QVV 



dJ„ 



< a z >-- 



dh 



(3) 



The non-analytic behavior and finite-size scaling of two- 
body NN entanglement of ITF model have been stud- 
ied intensively^. Then first, we look at the ITF case 
(J x = 1, J y = 0). In order to quantify the connection be- 
tween the non-analytic behavior of Entanglement when 
the system crosses the critical point and statistical me- 
chanics, we look at the derivative of the NN correlation 
functions and TS as a function of h (for simplicity we 
have consider the absolute value of the ground state en- 
ergy). 

In FigfT] we have shown the derivative of G yv with re- 
spect to the magnetic field versus h. 

As expected, there is no divergence for finite lattice 
sizes, but for the infinite chain dG vv /dh diverges on ap- 
proaching the critical point. A more detailed analy- 
sis shows that the position of the minimum QiMin) of 
dG vv / dh tends towards the critical point as hMin = h c — 
N~ e with 6 = 1.50 ± 0.03 (Fig.©, inset). Moreover, we 
have derived the scaling behavior of \dG yv / dh\h Min ver- 
sus N. This is plotted in Fig. ([2]), which shows the linear 
behavior of \dG vv / dh\h Min versus ln(N). The scaling be- 
havior is \dG vv /dh\ hMm = rln(A) with r = 0.25 ± 0.03. 
It is interesting that it is very close to the ITF correlation 
function exponent rj = 4. 

For an infinite system, dG yy / dh diverges on approach- 
ing the critical point as d( ^ h = /iln(|/i c — h\), where 
shown in Fig.©. Here fi is 0.27 ±0.03 but the analytical 
calculation shows that it is not completely independent 
ofh. 

According to the scaling ansatz 1 — , the two-site corre- 
lation functions, considered as a function of the system 
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FIG. 1: (Color online) Evolution of the — - 
different system sizes for J x = 1, Jy = 0. 



versus h for 
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FIG. 2: (Color online) Scaling of the minimum of ^ for 
systems of various sizes. Inset: Scaling of the position (hMin) 
of d( ^ vy for different-length chains. 



size and the coupling, is a function of N x ^ v (h — hMin)- 
In the case of logarithmic divergence which obtained 
in FigH it behaves as (dG yy /dh - dG yy /dh\, lMz J ~ 

F(N 1 /"(h-h M in))- 

In this way we have plotted 1 — exp(dG yy /dh — 
dG yy /dh\ hMa J versus N^ v (h - h Min ) in Fig.®. The 
curves which correspond to different system sizes clearly 
collapse on a single universal curve. Our result shows 
that v = 1 ± 0.001 is exactly correspond to the correla- 
tion length exponent of ITF (v = 1). 

A similar analysis can be carried on G xx , G zz — and TS 
(x z )- Our calculations show that the non-analytic and 
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scaling behavior of NN correlation functions and TS are 
the same as each other. The logarithmic divergence of TS 
near the critical point has been shown in inset of Fig.(|3]) 
in which /x = 0.27 ± 0.03 equals the scaling coefficient of 
dG yy /dh. In order to confirm our findings, we investigate 
analytically the derivative of correlation functions and 
TS using the exact ground state expression of ITF model 
which is obtained using the Eq.([2]) for J x — l,J y = 0. 
It is not hard to see that the ground state energy per 
site can be written as a complete elliptic integral of the 
second kindii. 
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In the vicinity 
— , and therefore 
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in the vicinity of critical point (A c = h c = 7 = 1) the 
ground state energy contain the nonanalytic contribution 
E G = (h-h c fhx\h-h c \. 

Similarly the G xx correlation function and transverse 
magnetization (< a z >), which are obtained using the 
derivative of ground state with respect to the J x and 
the magnetic field, could be written as a functions of 
complete elliptic integrals of first and second kindsi^, 
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^*(-,a) 

1 + A v 2' ' 



8 1 p( n < 



(5) 



where K(^,a) is complete elliptic integrals of first kind. 
Then after some complicated calculations, we could ob- 
tain ^L- ~ Snr ~ X z ~ (iki\h — h c \ in which the coeffi- 
cient of \i is function of h. 
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FIG. 3: (Color online) Logarithmic divergence of | ^ fa | close 
the critical point for the infinite system. Inset: Logarithmic 
divergence of transverse susceptibility (x z ) near the critical 
point for the infinite system. 
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FIG. 4: (Color online) Finite-size scaling of dG yy /dh for dif- 
ferent lattice sizes. The curves which correspond to different 
system sizes clearly collapse on a single curve. 
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dGpdg =\L 2 Ln(\h.-h\) 
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FIG. 5: (Color online) Derivative of next nearest neighbor 
correlation function (Gf 1 ) versus h for infinite size system. 
Right inset: Logarithmic divergence of next nearest neighbor 
correlation function near the critical point for the infinite size 
system. Left inset: Derivative of third nearest neighbor cor- 
relation function {G% x ) versus h for the infinite size system. 



III. SECOND AND THIRD NEAREST 
NEIGHBOR CORRELATION FUNCTIONS AND 
THEIR SCALINGS 



The numerical results shows that the first derivative 
of next nearest neighbor (NNN) correlation functions 
with respect to the magnetic field, diverge at the crit- 
ical point (Fig.®) while the second derivative of NNN 
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FIG. 6: (Color online) Evolution of the transverse susceptibil- 
ity versus h for different system sizes for J x = 0.75, J y = 0.25. 
Inset: Finite-size scaling of \ z f° r different lattice sizes. The 
curves corresponding to different lattice sizes collapse to a 
single graph. 



concurrence presents such non-analyticity 3 . This differ- 
ence could originate from the properties of the entangle- 
ment. The entanglement is linear function of the cor- 
relation functions and shows the collective behavior of 
them in which the details of the subsystem's behavior 
could be hidden, the logarithmic divergence of NNN 
correlation functions has been shown in right inset of 
Fig.© with fX2 — 0.51 ± 0.03. Despite the discrepancy, 
both of NNN concurrence and NNN correlation functions 
manifest the logarithmic divergence and same finite size 
scaling (y — 1). Moreover the third NN correlation func- 
tions's derivative diverge at the critical point (Fig®, left 
inset) and their finite size scaling behaviors are the same 
as NN and NNN correlation functions. It is important to 
mention that even when entanglement is not present (for 
ITF case entanglement is indeed completely absent for 
sites farther than second nearest-neighbors^), the QPT 
can be clearly revealed through the singular behavior of 
correlation functions. 

Second, we consider the case for J y ^ 0. Fig.® spec- 
ifies the divergence and finite size scaling of TS (Fig©, 
inset) for J x — 0.75, J y — 0.25 (correspond to the case 
7 = 0.5 in ref4). In this case scaling is fulfilled with the 
critical exponent v = 1 in agreement with the previous 
results and universality hypothesis. However the scal- 
ing of the maximum's position (huax) of TS has been 
shown in Fig.© which manifest the linear behavior of 
ln(h c — liMax) versus ln(iV). Inset of Fig.Q shows the 
scaling of the maximum of \ z versus the system size. 

We have also studied the correlation functions and TS 
of quantum compass models in the presence of homoge- 
nous and inhomogeneous magnetic field. This model 



FIG. 7: (Color online) Scaling of the position (htiax) of \ z 
for different-length chains. Inset:Scaling of the maximum of 
X z for different lattice size. 



shows a rich phase diagram which includes several criti- 
cal point depending on the exchange coupling 1 -. In this 
model the correlation functions of the various component 
of the two spins on odd bond are different from those one 
on even bond. Therefor this model is reasonable model to 
examine the above results. We have been able to obtain 
the divergence and scaling behaviors of the NN correla- 
tion functions and TS of this model 16 . This results justify 
the above finding. 

IV. SUMMARY AND DISCUSSION 

To summarize, we have investigate the properties of 
the correlation functions and we have shown that be- 
cause of their nice scaling properties, phase transition 
point can be determined from small systems with consid- 
erable accuracy without pre-assumed order parameters 
even for the cases where the pairwise entanglement is ab- 
sent. Then any smooth function of correlation functions 
(concurrence, quantum discord, reduced density matrix) 
could bear the same dominating long range physics about 
the phase transition as correlation functions do. How- 
ever, we show that the peak anomaly of NN concur- 
rence is exactly dominated by short range interactions 
of the Hamiltonian. So studying the divergence and scal- 
ing properties of correlation functions enable us to study 
the none-analytic and scaling properties of entanglement 
without direct calculation of entanglement (concurrence) 
which is elaborate calculation. However, for some cases 
in the one dimension, the correlation functions' predict- 
ing power of phase transition owe to their long-range be- 
haviors could be more reliable than same predicting by 
the entanglement where shows short-range behavior and 
divergence in its second order derivative. Moreover, our 
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results suggest that because of nice scaling properties of 
correlation functions, they are proper quantities to ob- 
tained the QCP from the finite size data where definition 
of a suitable order parameter (susceptibility) is difficult. 

Further investigations including blocks and multi-body 
correlations functions may be interesting to establish a 
precise comparison between universality and scaling be- 
haviors of correlation functions and entanglement at crit- 
ical points especially in higher dimension. Such topics are 
left for a future research. 
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